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COMPRESSIBLE EFFECTS AND MODELING 


Joseph G. Marvin 
Ames Research Center, NASA 
Moffett Field, California 94035 


t . INTRODUCTION 


The rapid development fif faster*, larger computers has been paralleled by 
an equally rapid development of computational aerodynamics. In light of this, 
it has been suggested that the computer and the wind tunrtel will perform equal 
and complementary roles in the development of new aerospace vehicles. To 
make that suggestion a reality it will be necessary that the understanding of 
turbulence be enhanced and that some practical means of modeling it in compu- 
tations be provided. 

Table I is a simplified summary of the status of computational aerodynam- 
ics that illustrates its rapid development, particularly since the start Of 
this decade. The stage of approximation of the governing equations has been 
divided into four progressively more complex categories, culminating with the 
viscous, time-dependent Navier-Stokes equations. Beginning in the 1930’s and 
progressing through the 1960's, inviscid linearized theory in various refined 
stages has been used in the design of many aircraft; however, limitations in 
this theory required that much of the configuration design be accomplished 
experimentally. In the 1970' s, development of inviscid nonlinear theory 
advanced more rapidly and is now nearing completion. 

Computations for transonic and hypersonic flight have been made for 
realistic aerospace vehicle geometries (notable examples are refs. 2-4). The 
major limitation of these computations is that they cannot handle separated 
flows. To provide that capability, the next stage of equation approximation 
requires utilization of the Navier-Stokes equations or their approximation. 

For turbulent flows, which are ubiquitous and very important in most aerospace 
vehicle applications, computations employing the time-averaged Navier-Stokes 
equations are in the early stages of development. The limitation of these 
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compute l ions is the accuracy ef the turbulence lik'd o 4 used to complete the sys- 
tem Of governing equations; this item now paces the development:. 

Once this stage of approximation reaches the point of practical utility, 
it is expected that the development of computations using the complete* Navier- 
Stokes equations in time-dependent form will begin. However, because the 
resolution scale for the smaller turbulent eddies precludes the use of any 
present day computers, an advanced computer is essential before practical com- 
putations can be made available. Thus, within the next decade, it would be 
possible to numerically simulate the flow about aircraft, including important 
Viscous effects, if the ability to model turbulence can be improved. 

The Navier-Stokes equations are the basic governing equations used to 
describe most fluid mechanics phenomena. They also apply to problems involv- 
ing turbulent flow where turbulence appears as a random, dissipative, three- 
dimensional phenomenon that involves many characteristic scales. To avoid the 
difficulty of describing every discrete turbulent motion possible* some type 
of averaging must be employed. Of the many possibilities, Reynolds time— 
averaging^ 5 ) has proved the most successful. Time averaging eliminates some 
information contained in the original equations, however, and also results in 
more unknowns than governing equations through the introduction of apparent 
Reynolds stresses for the actual transfer of momentum by velocity fluctuations. 
Deriving additional equations for those Reynolds stresses only results in the 
introduction of additional unknowns. 

Therefore, the extra unknowns must be represented by physically plausible 
combinations of quantities for which transport equations are expressed in 
terms of constants or empirical functions considered as known or expressible 
in terms of the mean variables. The problem of reducing the unknowns to equal 
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the number of equations Is referred to as the closure problem; the process of 
expressing the unknowns as transport equations in terms of empirical functions 
or Constants is referred to as "turbulence modeling." 

Historically, progress in turbulence modeling has been slow and deliber- 
ate, and has relied substantially on a few carefully controlled experiments 
performed over a range of test conditions. Such experiments could be called 
"building-bloGk experiments" because they provided the gage for establishing 
the credibility of computational techniques and even more importantly because 
they provided physically meaningful concepts that were used to enhance heuris- 
tic modeling ideas. The conference On the computation Of turbulent flows held 
at Stanford University in 1968< 6 ) used these key experiments to assess pro- 
gress in predicting incompressible, attached, turbulent flows. Later that 
same year, the conference on compressible turbulent flows held at Langley 
Research Center concluded that very few, if any, compressible flow experi- 
ments in the building-block category were being performed. 

If progress in modeling for compressible flows is to be made, it will 
also come through combining a broad experimental effort with developments in 
computational techniques and modeling ideas. For compressible flow, however, 
progress is even more difficult to achieve because compressibility introduces 
additional complexities, not only in the modeling concepts, but also in the 
experiments. The complete concept of turbulence modeling has three elements: 
experiments, intuitive modeling concepts, and the computer code (sde fig. 1). 
Each element is essential, and the elements are not easily separated* 

The process could start with the modeling concept or with the experiments. 
Historically, the process began with experimental observations that later led 
to modeling concepts. This trend is beginning to change in that model 


development and experiment: are being performed in parallel and coordinated 
efforts. Once the mode 1 ing concepts are established, the computer code can be 
assembled. This particular sequence is Essential for complicated equation 
systems because the modeling concepts can often alter the. ord£r of the equa- 
tion system or method of solution. Once the code has been established, it can 
be compared with and Verified by the experiments. If the experiments provide 
enough detail, they can guide changes in the modeling concepts and the process 
can be continued until the predictive capability of the code is established. 

An important aspect of the computation code development is that it be directed 
specifically to the geometry of the experiment and that it use exact experi- 
mental initial conditions so that no doubt can be cast on the comparative 
results . 

The status of modeling for compressible flows carl.be broadly summarized 
as shown in figure 2. For constant pressure flows, that is, flat plates or 
Uodies of revolution where the boundary layers are thin aiid the outer inviscid 
flow can be described, adequate turbulence models are available and there are 
apparently sufficient data to verify them. Pressure gradients introduce com- 
plexities that begin to test our ability to model turbulence* For mild 
gradients, adjusting, the constants in the models that handle constant pressure 
flows will usually suffice. For severe pressure gradients, however, models 
are not generally available and there is only a limited range of quality data 
to verify models now under development. Turbulence models ate not yet avail- 
able for the case when pressure gradients are strong enough to separate the 
flow, arid there are only a limited number of quality experiments for 


verification. 


In the course of this presentation, A Concept of turbulence modeling for 
Compressible flows will be reviewed. This review will provide an up-to-date 
status of the problem and an assessment Of the progress being made to solve it, 
A description of turbulence modeling will be given first, and various modeling 
concepts will be introduced. Because turbulence modeling relics substantially 
oh experiment, the next section will define what an adequate experiment is, 
and what the present capabilities of measurement techniques are. In the final 
section, the Status of turbulence modeling for two-dimensional flows will be 
discussed and examples will be giver, for attached and separated flow problems. 


II. MODELING CONCEPTS FOR COMPRESSIBLE FLOWS 
Before describing the various concepts for the turbulence modeling used 
later in describing the status of modeling for compressible flows it is worth 
while to show briefly how they were developed from the basic governing 
equations* 


A. Governing Equations in Mass-Averaged Variables 
The Navier-Stokes equations for a compressible, heat conducting, perfect 
gas may be written as follows: 
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Where* the stress tensor t j j arid the heat flux vector q j are given by 



the bulk vlscositv X is equal to -2/3 p, where p Is the dynamic viscosity; 
the thermal conductivity is k; . is the Kroneckfir delta, equaling 1 when 
i = j and equaling 0 When i ^ j . A summation is implied when indices are 
repeated. 

As is well known, when dealing with flows that are turbulent, the solu- 
tion to equations (1) through (3) becomes practical only if one employs some 
type of averaging. Following Reynolds, for example, the equations are aver- 
aged over a time that is long compared to the time associated with the predom- 
inant frequencies of the turbulence. For compressible flows, however, mass- 
weighted averaging is often used, father than the usual time averaging, 
because the compressible form of the equations can be reduced to a form analo- 
gous to the incompressible form (see, e.g., ref. 8). 

the advantage of mass averaging can be illustrated as follows. In the 
usual averaging procedure, velocity is usually written 

u i (x i ,t) - u i (x i ) + u^(x 1 ,t) (6) 


where the bar denotes the time average of u'^Cx^t). the time average of a 
quantity is obtained from 
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Introducing the concept of time averaging into the continuity end momentum 
equations, the following equations result: 
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(f . u jP -u” - pu«uj - P"uju«) (9) 


Note that for incompressible flows dp = 0 and the equations reduce to the 
familiar incompressible form. However, aside from the usual problem of sup- 
plying additional relationships for the mean effects of turbulence, for 
example, pup 77 , the compressible form of the equations contains fluctuating 
density-generated terms such as p"uj\ This complexity can be formally avoided 
by introducing mass-weighted averaging. 

Define the following mass-weighted velocity 

u* - pTL/p 


( 10 ) 


where the bar denotes the usual time-averaged value and the tilde denotes 
mass-weighted averaging. The instantaneous velocity can be written 

Ui<Xi,t) ■ ^i( x i) + u J/ x i»t) t 11 ) 

where u| is the Superimposed velocity fluctuation. The relation between 
u^ and u^ is obtained as follows: 

pujtx^j - (p + p'^u^Cx^ = 0 

or uj(x ± ) = -p"u[(x i )/p 

Taking the average of equation (11) and rearranging 
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follows that 


u j (x i ) - u i (X i ) = p"u"/p 

Thus the difference between the mass-weighted and time-aVeraged Velocities 
depends 6h a density-velocity correlation. Similar relationships exist for 
tettperature and enthalpy. Later, relationships between various shear stress 
and heat flux quantities in mass-weighted and time-averaged systems will be 
needed. Without derivation, these are written as follows: 
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Substituting mass-averaged variables into equations (1) to (3) yields: 
Continuity 
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Equations (3.3) to (15) now have the same form as the laminar Navi er-Stokos equa- 


tions except that the averaged variables appear throughout and the so-called 


turbulertt Reynolds stresses and heat fluxes appear. The equations are simpler 

than those given when the usual time averaging Is introduced. These equation!? 

alsd have the same form as the incompressible time-averaged equations except 

that the so-called Revnolds stresses, -ou!u|, include fluctuations in density. 

i j 

This latter fact does indeed Complicate the process of modeling the turbulence. 

The equations in their averaged form have additional unknowns: the tur- 

bulent Reynolds stresses and heat fluxes. The problem of reducing the number 
of unknowns to equal the number of equations is called the closure problem. 

The process of expressing the unknowns as transport equations or functions in 
terms of known quantities is called turbulence modeling. 

The governing equations can be supplemented by additional equations for 
the turbulent kinetic energy and the various Reynolds stresses. ^ Following 
the development in chapter 2 of reference 9, these may be written as: 

Turbulent Kinetic Energy 
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Inspection of the conservation equations shows that the number of unknowns 
exceeds the number of equations. Thus, the next step in arriving at a closed 
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Hot of o qua 6:1 on 8 will bo to tiiodol certain of the tefros. Tills ran. bo arcom- 
p I is l led in progressively more complex stops in the sanio manner as for Jnc.om 
p ro s s I h .1 o f X own ♦ 


B. . Modeling Concepts 
I . An Alge br aic Clos ure Model 

A straightforward technique that has received much attention and verifi- 
cation for a broad range of boundary-layer flows is to express the turbulent 
shear stress and heat flux in terms of the mean flo.w variables. .Bousinesq ' s 
eddy viscosity concept is introduced. Define the turbulent Reynolds stioss 
tensor as 


then let 





( 18 ) 


where t ..is considered a scalar quantity referred to as the eddy viscosity 
because equation (18) now has the same form as the laminar stress tensor, x ^ . . 
The problem of modeling is now reduced to one of describing the behavior of 
the eddy viscosity in terms of the mean flow variables. Several techniques 
for accomplishing this when the flow_is .compressible are described below. 
Essentially they are extensions of techniques used for incompressible flows. 

A successful tflodel was developed by Ccbeci and Smith. It uses Prandtl s 
mixing-length hypothesis (i.e., t = H 2 [ O&j/axj)- + O&j/axp]) in the inner 
region with Clauser's eddy viscosity in the outer region tc arrive at a two- 
layer model. Also, in order to appropriately model the region of the sub- 
layer and the inner log-law region, vrift Driest' s damping function is 
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introduced. The model has been successfully applied to many compressible 
equilibrium turbulent boundary-layer flows. 

Inner region 

( Du. auA 

5IJ + * h ) ( 

Outer region 

i = 0.0168 p u 6 */l + 5 . 5 (y /6) 6 (20) 

outer e e k 

where 

D = [1 - exp(- y/A)] (21) 

and 

A * A + v(t w /p‘ t ) (22) 

and 

/•^edge 

“ J (1 " u/u*)dy (23) 

Following Morkovin's hypothesis that compressibility does not affect the 
description of the length scales, the constants in equations (19) to (23) are 
given values derived for use in incompressible flows. In equation (22), A + 
contains a reference density which the originators assume to be the local den- 
sity across the sublayer and in most applications the wall density is used to 
represent this reference value. The value of A + can be adjusted to account 
for pressure gradient and mass transfer at a wall in the same way as for 
incompressible flows. (See chapter 6 Of ref. 9 for a discussion and suitable 
derivation of eq. (22) tOr compressible flows.) Hereafter, .this model with 
A + » 26 will be referred _to as the. baseline, turbulence model. 
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For £l.d>wfi with compressibility 16 is also required that the- ‘‘rtergy equa- 
tion be solved . To express the turbulent he.1t flux ah eddy eonduetivity is 
introduced and a turbulent Praridtl number is formed 

Pr t - f 

e h 

Expressions for the turbulent Prandtl number have been developed for boundary- 
layer applications. (See, e.g., chapter 6 of ref. 9.) 

2, An Algebraic Closure Model with Relaxatii^n 

The baseline turbulence model is not always successful. For equilibrium 
turbulent flows with mild wall pressure gradient, modifications to A + will 
usually result in adequate predictions, but for flows with rather Severe pres- 
sure gradients, such as those ahead of separation where the flow is not in 
turbulent equilibrium, modifications to the model are required. Recently, 
several studies have shown limited improvements for these situations can be 
obtained by introducing the concept of relaxation* See, for example, refer- 
ences 10 and 1.1. An additional differential equation for eddy viscosity has 

been written as 

de = C ~ Ee q (25) 

dx L 

where L is some relaxation length Snd e eq is the eddy viscosity given by 
the usual prescription, for example, equations (19) to (23). Some authors 
have applied the relaxation concept to the outer eddy viscosity only, arguing 
that the smaller eddies in an inner region should be able to adjust immediately 
to Changes in the mean flow/ 10 ’ 11 * Another author applied the concept across 
the entire boundary layer/ 12 ) In either case, the relaxation concept is 
introduced to account for the experimentally observed fact that turbulence 


Cannot adjust instantaneously to changes in the mean flow. However, this 
particular relaxation concept has not beea successfully generalized and the 
relaxation length scale, L, usually expressed in terms Of boundary-layer 
thicknesses, is not universal. Examples employing this Concept will be shown 
later. Hereafter this model Will be referred to as the relaxation mode).. 


3. A Single Transport Equation Model 

A more straightforward approach to modeling for nondquilibrium turbulent 
flows is to introduce the transport equations for Reynolds stresses. The 
approach is much more complicated than for incompressible flow, however, 
because terms involving compressibility must also be modeled. Little experi- 
mental evidence is available to help in this latter process. To illustrate 
the problem, the development of these additional modeled terms for the ohe- 
equation model derived in reference 13 is briefly described. The eddy vis- 
cosity in this application is formed from the product of an algebraic length 
scale, expressed in terms of a turbulent Reynolds number, and a velocity 
formed from a solution of the turbulent kinetic energy equation. The method 
is an extension of Glushko's incompressible model. 

A more general scalar eddy-viscosity relationship, described by Hinze/ 15 ^ 
wherein the average turbulence pressure is extracted from the stresses, can be 


written for the turbulent stresses. 
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rm rotlucing equation (2ftn) into oqu.it. I on (lb) and rest ruotur I ng tin' Lt-rnis 
involving vo .lot: i t y •* p r e:» 8 u to gradient 1 , and veloei ty-slu-af stress correlations, 
tl.e following tnrbuloiit kinetic; Ortorgy equation is obtained: 
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The turbulent kinetic energy equation in the mass-weighted averaged system has 
the same form as the incompressible equation but the problem of modeling the 
terms on the right-hand side of equation (27) differs from the incompressible 
approach because extra pressure velocity correlations and molecular-shear and 
velocity correlations are introduced through the mass-weighted averaging 
process. (See ref. 13 for a complete discussion.) To exemplify this, con- 
sider the pressure-velocity correlations. Substituting p P + P into 
those terms results in the following group of terms 


9u! 3u^ 

W ‘ - Tx7 (u i p) * u i p ” + J »h + p ” 


(28) 


Both "uT and its derivative are zero for an incompressible flow and only the 
second term on the right-hand side remains; it is usually grouped with the 
second term on the right-hand side of equation (27) and modeled in terms of a 
gradient of k. (See chapter 5 of ref. 9.) The choice for compressible flow 
is either to assume that these remaining correlations are negligible or to 
model them in such a way that they tend to Zero when the flow is incompres- 
Fol lowing the procedure of reference 13, rewrite equation (28) as 
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14 


sible. 


Equations for have been obtained for two cases; Otte based oft the assump- 

tion of constant total temperature within an eddy and the other based on the 
assumption that the state variables in an eddy behave in a poly tropic fashion. 
In either case, generality is maintained by assuming that pressure fluctuations 
can be expressed in terms of density and temperature fluctuations. Assuming 
that the gas behaves in a polytropic manner, then 
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where n is the polytropic coefficient and terms Of second Order are 
neglected. Using this expression to relate pressure and temperature, Rubesin 
arrived at the following equations 
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where a is a speed of sound and £ is a coefficient also containing the 
polytropic exponent. Thus, it is possible to arrive at expressions for the 
additional correlation terms introduced by compressibility , but they must also 
be modeled through the introduction of values for n and £. The appropriate- 
ness of the values suggested by Rubesin, however, has not yet been demon- 
strated. Experimental measurements of the terms are not available* and the 
model is only now being tested extensively against experiment. 

The complete equations describing the one— equation model are given by 


equations (75) through (91) of reference 13; they will not be repeated here 
Examples shown later refer to this as the l-equati_o_n_jriod el. 
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4. A Two Equ at ion Model . 

A more generalised modal can be developed by loaning an eddy viscosity 
from a velocity and length scale that ar& both described by transport equa- 
tions. For example, the previous one-equation mrtdtil could be generalized to a 
two-equation model by providing one additional transport equation for the 

length scale as originally proposed by Rotta/ 16 ^ 

Another approach that is receiving close attention now and for which some 
results will be presented later is the two-equation model recently developed 
by Wilcox^ 7 ) for the boundary- layer approximations to the Navier-Stokes equa- 
tions. This model uses an energy equation and a dissipation-rate equation 
that is a modified form of Saffman’s pseudo-vorticity equation. Compressibil- 
ity has been introduced through mass-weighted averaging as described pre- 
viously. The modeled form of the equations developed in reference 17 are 


\ (pk) + IjxT (pa j k) = fy “ 3 *“) k + ¥ 


(y + o*pe) 


9k 

9y 


J - C*pk 


dUy 

9x, 


(32) 


and 


*[• W * b (?)'] ■ I 


- 9u 

ap 97 


6 + 2a 




+ 


3y 


(y + ape) 



(33) 


where the term containing C* is an additional term that includes the dilata- 
tion of the mean flow due to compressibility and where the eddy diffusivity 

(i/p) is 

- p£ 


In their latest interpretation, Wilcox and Tracy identify k with 
9/4 V'~ 2 , rather than the turbulent kinetic energy. Thus, inherent in their 
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model interpretation is the arsumpt. ion that ( <u' ) > ' • < <w' >)'■ ; ( <v’ > ) Z =* /|:3:2< 

The examples to he presented ’.iter vl-U handle ine add! t tuna 1 'Oi;i|ies:;lbl< 
modeling terms K and tn by choos ! hr zCivo for both Ond thus reducing tin- equa- 
tions to their form that is usually associated with ineoiiipresiulile llow. Ihus, 
compressibility is accounted for only through introduction oi i he mean density. 
Examples given later refe r to this model as the Wj i cox-Tracy . 


5 . A Complete Reynolds Stress Model 

Applications oi the modeled form of the *'.,11 Reynolds stress equations 
have also been attempted for compressible flows. See for example reference 18 


where Donaldson’s invariant modeling ideas ate applied to a compressible >1o.’. 
Donaldson's approach has been receiving focused attention recently i ct 
boundary-layer flows and some examples are presented .latei . Rather than 
employing mass-weighted Variables, this approach begins with the complete 
time-averaged Navier-Stokes equations and closure is accomplished by modeling 
the third- and fourth-order correlations in terms of the mean variables and 
second-order correlations according to the prin iples of invariant modeling. 
For the sake of brevity, the equations and the modeled terms will not be given 
since they Comprise a set of 12 equations and .1.2 unknowns, the reader \u 


referred to references 19 and 2u i.:>c a 
the equations. For the examples Lo bo. 
wore taken from reference 19 and were 


discussion and formal representation oi 
given later, all the modeling constants 
hot adjusted. Later , this mode l i 


referred to a the ARA P no.lcylj 


ITT. EXPERIMENTS FOR TURBULENCE M'''DE' ’’"T 
Turbulence modeling is heuristic by nat" e there fore t 11 : o:.t ■ 

give pxpeHnvnital data for guidance and verify tion. V He 'h re is a rat 1 


17 


T 


extensive <lnin 1 t o support modbllng for i ncompress 1 T> l e f ! nv/s; j the duLa 
base for .support i ng compress lb l.e flow modeling, Is ralbe.r small. Only ■•ec-.-nily 
I, as chore boon a concerted attempt to rectify this imbalance. < 2 1 ^ A s i )‘,n i I i •• 
cant contributor to this rather unfortunate imbalance is the diflieulty .>) 
measurement when Mach number must be considered. Jh this section the elements 
of experiments to support modeling will be discussed, the status oi. measure- 
ment techniques will be described, and an experimental foundation for develop- 
ing a model, that will be useful for computing two-dimensional flow fields will 
be proposed. 


A. Elements of Experiments Supporting Modeling 
Figure 3 shows the Mach-Reynolds number domain for aerospace vehicles and 
it gives an indication of the range of conditions over which adequate turbu- 
lence modeling must be provided if the concept of numerical simulations of 
flow fields over complete aircraft is to become a reality. The upper limit on 
Reynolds number based on vehicle length ^ 22 ^ is shown as well as the upper 
limit based on a typical wing chord, Mach number varies between subsonic and 
hypersonic, encompassing the range encountered by commercial passenger 
Vehicles and NASA's space shuttle vehicle. Experiments that verify and guide 
the modeling process will have to be performed over this range of conditions. 

Experiments that support the modeling process can be classified according 
to the type of closure proposed. Bradshaw^ 2 ^' broadly classified these 
closures as first and higher- order . First-order closures correspond to 
closures where second-order correlations like the Reynolds stresses are 
expressed in terms of first-order correlations like the mean velocity, namely, 
algebraic mixing length or eddy viscosity formulations. Higher-order closures 
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refer Co those in which third- and lUgher-order r.omlat inns are expressed »<•< 
terms of second- and higher order correlations l*.y Inf.foduc i ng adtlh nmal 
appropriate transpor; equations . 

Table II gives the elements of the experiments required, depend i ng (, n the 
type of equation closure. Bv their nature, verification experiments require . 
documentation of mean and surface quantities over the practical ranges of 
flight Math and Reynolds numbers. Those experiments are useful for any clo- 
sure technique used. First-order mod- ling experiments arc those that require 
measurements of the shear stress and heat fluxes across the flow Held because 
these quantities provide insight into concepts used to model these terras and 
provide closure. 

Such experiments can, but need not be, attempted over the. -complete Mach 
and Reynolds number ranges because the verification experiments can test the 
ability of the model tc perform outside the domain where, for example, the 
shear stresses have been measured. But they must be performed at Reynolds 
numbers sufficiently high to establish fully developed turbulence. Higher- 
order modeling experiments require that fluctuating measurements be made across 
the flow field. Depending on the order of the closure, more and more informa- 
tion on the fluctuating field must be ascertained, as was the case for fi’-st- 
brder modeling experiments, data need not be acquired over the complete Mach 
and Reynolds number domain, provided that verification experiments are avail- 
able. Tdeallv , one well-conceived experiment could suffice for all three 
types listed in table II. Mo cover, the same <.-p.-riment.- l apparatus and 


instrumentation can be used to eliminate experimental uncertainties. Colo 
emphasised this latter aspect- v.n<-n commentin' oh the 1 " t '" ; USCl1 tVu ,Jd ‘ 
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measurement, complete documentation for future fbfetertoc, and the eliroi. nation 
of throO-dimCnsiona.1 uncertainties by, for example, testing axi symmetric 
configurations. 

As mentioned previously, there are few experiments suited to guide and 
verify turbulence modeling for compressible flows. Figure 4 compares the 
domain of available experiments that can be used to test modeling concepts 
with that for vehicles. The shaded area represents the domain where partially 
documented experiments were performed for flows where zero- or mild-pressure 
gradients were impressed on the flow. Reference 25 presents a table summariz- 
ing most of these experiments and their measurements. In almost all cases, 
the experiments were performed at low Reynolds numbers and their decu..ientaLxon 
was incomplete; that is, some experiments measured only a single quantity and 
no single experiment was complete, at least not. insofar as ir, admg all the 
elements listed in table II is concerned.. It has been d'rficult for analysts 
to use these experiments to predict confidently the trends in Reynolds numbers 
for actual flight conditions; nevertheless, there appear to be sufficient data 
of the verification variety to make definite conclusions regarding the ability 
to model turbulence for attached flow boundary layers neat equilibrium as will 
be shown subsequently. 

The unshaded area in figure 4 represents the domain of more recent experi- 
ments where shock- induced separation was studied and where specific, attempts 
to document the complete flow field features, including the turbulent fluctua- 
tions and their correlations, have been made or are under way. These latter 
experiments and companion computations, which together are being used to 
establish adequate turbulence models, are discussed later. 
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Before discussing ;iriy of t-hcsA experiment s * L'lio question of what roll he 
accompli shed experimentally in compress I bio flows where the dl fflrnlt c • • • : i 
enviroiimeiiC of high Reynolds number.'! and high Mach numbers is pivsem tin i * t be 
addf cssf'd . Over the last several years, a considerable amount <T i nstrum n< 
development lias been undertaken that tends to make c:n •■ctatidns optimist 1 
in this regard. 


B. Status of Measurement Techniques 


1. Mean Flow Measurements 

Meaa flow measurements of Hie type required for verification experiment:; 
will not be discussed in detail. It suffices to Say that in most instance: 
means are at hand .to measure mean-flow quantities Such as velocity, pressure 
and temperature, and wall values of heat transfer It is worth mentioning, 
however, a recent development in the technique for measuring skin '-'ricrinn 
which heretofore has been difficult in cases where pressure gradients are 
applied. 

Ludweig's concept of measuring surface shear with a heated surface ele- 


ment has been reintroduced recently by RubetHn et al. 


lor use m compres- 


sible flows With rather sever., pressure gradients. By imbedding small diam- 
eter wires (1 mil) in a low conductivity material contoured to the shape of 
the model surface, it has been possible to obtain reasonably accurate measure- 
ments of surface shear. These gages have the distinct advantage of being 
small in size, insensitive to pressure gradient, and they can be used whether 
the flow is laminar or turbulent. A disadvant go is that they vit be cal'’-, 
brated over a ran, a- of apnropr i c condition n situ, but o lar this has not 
been a major problem. Exan-p] os of h measurements will be given later. 




2 , Shfrar Si tosh Mo astir emcrits 

In addition to mean-flow measuremefit a shear stress distributions will be 
required for iirr.t-order modeling experiments . For xe.ro pressure gradient 
boundary- layer flows at supersonic speeds, moment um-b a lane e techniques using 
mean flow-field data have been used successfully to extract the magnitude and 
distribution Of the shear Stress. Using this technique, Sandborn was able to 
show that the appropriately normalized shear stress was essentially indepen- 
dent of Mach number and Reynolds number in the supersonic speed regime to 
M« 7 . ( 27 ) The results ate shown in figure 5 taken from reference 27. The 


shaded region represents the distributions from_lA different experiments 
covering Mach numbers between 2 and 7 and Reynolds numbers, based- on momentum 
thickness, between 2,500 and 12,000. At higher Mach numbers the picture is 
not so clear, but this may be a direct result of the fact that most of the 
higher Mach number data were taken on nozzle walls where the boundary layer 
had undergone distortion by severe pressure gradients in the throat regions. 
Any proposed turbulence model for zero pressure gradient flows should be able 
to predict the normalized shear distribution in figure 5. Another important 
contribution of this, universal shear distribution is that it can also be used 
as a standard against which other direct measurement techniques can he tested 
as will be shown subsequently. 


Unfortunately, when pressure gradients are imposed on the boundary layers, 
momentum-balance techniques that employ mean profile measurements ate usually 
not satisfactory for evaluating the magnitude of the shear stresses. This is 
particularly true when gradients that load to separation are imposed because 
the shear stresses comprise only a small part of the overall momentum balance; 
unfortunately, the shear stresses are the important stresses that determine 


the eventual State of tlie flow field. Therefore, methods for measuring the 
shear stresses directly must be developed. Hot-wire anemometers and laser 
veldeimetets have recently been used for this purpose and it is important to 
btl'fly review their applications and limitations. 

Two major problems arise when employing hot-wire anemometers in compres- 
sible flows where high Reynolds numbers and Mach numbers are present, the 
first of these is the practical problem of wire durability. The second and 
obviously most important is the measurement interpretation. 

Mikulla and Horstman^ 28 ) and subsequently Mateet artd Brosh^ 29 -* developed 
"backed" wire probes that successfully eliminated the problems associated with, 
wire durability. Sketches of the two probes taken from .these references are 
shown in figures 6 and 7. The dual-wire wedge probe and the separate cross- 
wire wedge probe were used to obtain measurements of mass flow and vertical 
velocity fluctuation, and their correlation. The triple-wire probe was used 
to obtain the correlation of total temperature and vertical velocity fluctua- 
tions. The dual-wire probes were operated with constant temperature anemometer 
systems; the normal wire of the triple-wire probe was operated with a constant 
current system. Frequency response was ...considered adequate (50 kHz)~. 

Backing the wires in this manner eliminated wire breakage and Strain 
gaging problems; it permitted operation of cross-wires at high overheats with- 
out causing sag and maintained the wires at 45° angles so that vertical veloc- 
ity fluctuations cduld readily be obtained when the probe was positioned in a 
plane normal to the wall. The main disadvantage of these backed wires was the 
possibility of heat conduction between the wires through the backing material; 
such conduction could affect the mass-flow sensitivity calibration. For the 
high Mach number applications of Mikulla and Horstman , (28) conduction between 
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i.lU' wires was »uf f i c I .mil to preclude a sensitivity calibration during siendy- 
3 Unto opera t ion because the wires were operated at Uenipei itui cs significant I./ 
Iiij'.hdf than the test total temperature. (1500°R). Mtknlla and llorstrnan wen-, 
however, able to inter the neiifSit i v:l 1 y cal i brat. I an by ope rat: ini’ the backed 
wire ill a known boundary-layer Clow previously measured with an unbacked normal 
wife, by relating, their previously measured mass- 1 low fluctuations to l lie 
voltage output from the dual wire, they could infer the sensitivity eoefi i- 
cients. Fortunately the probe of Mateer and Brosh^ 9 ^ was. operated at a lower 
temperature and was not affected by this problem; consequently, they were able 
to calibrate directly. 

Extracting turbulence information from hot wires operated in compressible 
flows can be accomplished in. a straightforward manner. However, certain 
assumptions are* inherent in the process and they require some explanation. 
Following Morkovin, ^ 30 ^ the expression for wire voltage fluctuations is written 
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S (pu) 


(pu)" 
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(v") 
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(T t a ) 


(34) 


where S ( ) represents the wire sensitivity to mass flux (pu) , wire orienta- 
tion angle (.’>), and total temperature (T t ) an i the double prime refers to a 
time varying quantity. The plus or minus sign ahead of the second term is 
determined by the sign. convention chosen for the vertical velocity. Morkovin 
showed that, equation (34) was restricted to wire operation such that M u> sin v 
was greater than 1.2, a range where the sensitivity of wires to density and 
velocity changes were equivalent and could be coMilvned into a single sinsitiv 

S (^)* 

More recently, a study (3l) similar to Morkovin* s but including more wire 
calibration data, shoved that equation (34) is also valid at lower Mach 


numbers provided the wire overheat, [t Wr - (f w - T f )/t r ], and wire Reynolds 
number, Re c , were above 0.5 and 50, respectively. The conclusions Of the 
study were based on calculations of the sensitivities obtained from the corre- 
lations of wife data by Berhdns^ 32 ^ and corroborated by comparisons with sen- 
sitivity data from normal wires obtained by Rose and McDaid.^ 33 ^ bxamples of 
thd results, given in figures 8 and 9 , substantiate the Conclusion that equa- 
tion (34) is also valid at lower Mach numbers provided the restrictions men- 
tioned above are met. 

Modal analysis techniques similar to those proposed by Morkovin^ 30 ^ for 
a normal wire can be used to solve equation (34) for the fluctuating quanti- 
ties because the sensitivities are known functions of wire temperature that 
are Obtained by suitable calibrations. Squaring equation (34) and time aver- 
aging yields a set of equations with six unknowns, the three fluctuating time- 
averaged variables and their cross products. Conceptually, a solution is 
possible if the wire is operated at a minimum of nix overheats. However, to 
obtain accurate values of the six unknowns extremely accurate values of the 
sensitivity coefficients and fluctuating voltages must be obtained. To 
increase the accuracy of the process a normal wire can oe '-sed to obtain 
(pu)" and T t " and their correlation independently an'.-, the yawed wire 
results only used to obtain the remaining unknowns. Rose used this technique 
to sort out the fluctuating quantities and their correlations in a turbulent 
supersonic boundary layer in the presence of an adverse pressure gradient. (34/ 
In Rose's case, this process was made easier since the terms involving 
and its correlation with v" .were of smaller order than the remaining terms. 
An alternate technique can be employed that does not depend on wire operation 
at a large number of overheats. For obtaining information on the fluctuating 
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m.i:is flow and vortical velocity and the it correlation, the wires may be oper- 
ated at high enough temperatures ho that they are sensitive mainly to mass- 

flow fluctuations, S/~, » S,7 p - V Conversely, they may be operated at low 

Ipu; u i> 

temperatures where they a if 6 Sensitive mainly to temperature fluctuations, 

s — >> s / — . (The latter application is difficult to achieve in practice 

(T t ) (pu) 

with "backed Wires" operated at constant temperature so they are. usually 
operated at constant current.) 

Consideration, must now be given to how shear stresses might be obtained 
with hot wires. Provided S^— ^ >> S > a dding and subtracting the signals 
from a dual-backed wire and time averaging will provide data on the root-mean 
square -values, <(pu)"> and <(v)">, and the correlation R^y.y* The 
measurements then carl be used to give 


(pu) "v" = R( p -- y, f ^ 7[ < (pu)" > 1 [ < (V") > 1 


(35) 


Usually an independent measurement of <(pu)"> is obtained by operating a 
single normal backed wire in the same flow. 

Expanding the left-hand side of equation (35), 

(pu)"v" = “u"v" + Op"v" (36) 

Introducing the polytropic gas law and neglecting higher-order correlations 
such as ? r f"/pf compared to ^/p and T^/T, the following equation ihay be 

written 
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(37) 


p - P n - 1 T 
Equation (37) can be written in terms of total temperature by introducing the 
perfect gas energy equation. 


o"v" h 1 { * "I'.Tt'i _ _y_ 1 
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(38) 
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Substituting equation (38) into equation (36) results in 


pu'V 


"Ip* 




.. (pu) 'V' - f Tt'V'j 


( 39 ) 


r r -[u 2 /(n - 1)]) 

Recalling that the shear stress of interest in the governing equations in 
section II was written in mass-averaged variables, the additional relation- 
ship is required: 


pu' v* 


. P'V’ P"v" 
P 


pu v 


(40) 


In many instances the second term on thd right-hand side of equation (40) can 
be assumed negligible compared to the first term. Equations (39) and (40) 


show how shear stresses can be obtained from hot-wire measurements of (pu)'V' 
and T^V\ The latter correlation may be obtained by operating a normal wire 
at low overheat such that » S ( — , with a constant current system 

along with a yawed probe operated at high overheat where S ( — ) » s ( f^)‘ 
the triple probe (26) mentioned previously was operated in this manner. Addi- 
tional manipulation of equation (39) to eliminate the temperature-velocity 
correlation in favor of the introduction of a turbulent Prandtl number can 
also be accomplished. Following the derivation in reference 35 the following 

is obtained. 

— rr-M _ fi + 1 mz\i - c _^/ - T 1 r(pu) 1 ^ (41) 

pu v - |1 + n „ 1 Pr t M ^ c p u(du/dy)JJ 

Expression (41) has been employed for both adiabatic and nonadiabatic 
wall temperature flows under the assumption of isobaric turbulence (n - 0) and 
for Pr t =0.9. Use of this expression does not require a measurement of 
temperature-velocity correlations. In reference 35, a comparison of pu’V 
with similar data using measured temperature velocity correlations showed 
excellent agreement. In all applications of hot-wire measurements to obtain 
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shear stress t he polyt ropi c. exponent h must he assumed. Pof t unat >• I > , 

Ktst let's measurements of adiabatic, boundary layers tu Mach number 5 showed 
that the IsoburiC relationship Cor turbulent eddies is adequate and so n *= 0 
is appropriate. No information on the exponent, n, is available above this 
Mach number, however. 

In principle, the laser velocimeter does not depend on the assumptions 
relating density and velocity fluctuations to temperature and velocity fluc- 
tuations. By intersecting two coherent light beams in a small volume within a 
flow field a precisely known interference fringe pattern is established in a 
plane perpendicular to their bisection. Any particle passing across the pat- 
tern alternately emits light .that can be picked up with a photo sensor. Given 
the known distance between fringes and the tittle duration of the light pulse, 
the velocity of the particle may be measured. Provided the particle dynamics 
within the flow field are understood, a flow velocity may be obtained directly 
with this instrument. The technique also has the advantage of being nonintru- 
sive. The problems of measurement interpretation should not be overlooked, 
however. Particle seeding techniques and measurement biasing toward higher 
velocities are but two of many problems encountered in the application of this 
instrument. It is not within the scope of this paper to discuss these prob- 
lems so it suffices to mention that problems do exist and that they are being 
studied. However, recent applications where alternate means of measurement 
are available suggest these problems are not significantly influencing the 
results in many applications. See for example, reference 36 where favorable 
comparisons between hot-wire measurements arid laser velocimeter measurements 


are shown. 


Turbulent shear stress may be obtained in a manner analogous to that for 
slanted hot wites.^ 27 ) The optical axis of the velocimeter is rotated so that 
the fringes of one pair of crossing beams ate oriented at a +45° angle to the 
flow direction. The velocimeter will measure a velocity (assume v = 0) 


V, - t(u + u") - v"] 

n 

Rotating the axis to -45°, the velocity measured is 

v ? = (u + u" + v") •“ 

/2 

The mean values of v^ and ar6 
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(41b) 


Where N and M 
variances become 
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are the total number of discrete velocity measurements. 


(42) 

The 


E 


i-1 


li 


- V 2 


N - 1 


and 


M 

E. 


(v 


2J 


- v 2 ) 2 


M - 1 


Squaring and taking the difference of equations (41) , multiplying by p 
results in 

5 u'V * (a 2 2 - a^) 


(43) 


(44) 


Recall that 


= £ TTV 7 + p'VV 7 

and therefore shear measurements with the laser velocimeter neglects the triple 
correlation. The shear stress in the mass-averaged variables must be obtained 


from equation (40) . 



Figure 10 shows the results from several experiments where shear stress 
measurements wfere made across turbulent boundary layers without pressure 
gradients* The boundary-1 ayo.r-cdgc Mach number In these experiments ranged 
from subsonic through hypersonic. The solid line indicates the expected 
variation of the normalized total shear stress that was shown in figure 5 and 
which is independent of edge Mach number. The measurements, obtained with 
hot-wire and hot-film anemometers and a laser velocimeter , show the variation 
of the turbulent component of Shear stress. Except for the decreases neai the 
wall (y/6 < 0.3), due mainly to the relatively large probe size and not the 
difference between the total and turbulent magnitudes of shear nor the neglect 
of triple correlations, the data agree reasonably well with the expected 
trends in shear distribution. Work is in progress to resolve the differences 
near the wall and some success had been achieved. For example, use of a 
split-film anemometer minimizes the probe scale effects and suppresses the 
decrease to a location much nearer the wall. Fram„s.uch data, mixing 

lengths or eddy viscosities .in the outer regions of a turbulent boundary layer 
can be determined for Use In first-order modeling concepts. 

3. Fluctuating Flow-Field Measurements 

Higher order modeling experiments require measurements of the individual 
fluctuating flow variables such as <u">, <v">, <w">, and <T > and perhaps 
their spatial .derivatives. As discussed in the previous section, velocity 
components may be measured directly by the laser velocimeter and under certain 
restrictions by the hot-wire anemometer. The restrictions on the latter device 
involve the choice of polytropic coefficient, n, to obtain streamwise velocity 
fluctuations <u" > and not on <v" > or <w" > , which may be obtained directly 

from a backed wedge probe if S^--^ >> ** (T«- ) . 


Choice of the isobarlfc value., n ■ o, involved In most Instance!? to 
Obtain <u">, U substantiated at M = 3,( 3f >) and can be trtfeffcd from 
Kisder's measurements to M - 5 . Temperature fluctuations may be obtained 
with the hot-wire anemometer. An example of fluctuating flow data obtairteu in 
a high subsonic boundary layer is .shown in figure 11. (these data were 
obtained with hot-wire and backed hot-wire anemometers.) The Solid lines are 
the usual ratios of 4:3:2 for ( <u' >) 2 : ( <w* >) 2 : ( <v‘ » 2 in incompressible 
flows. The measurements of the individual velocity components agree with the 
incompressible data. Density fluctuations at these relatively low Mach num- 
bers are small compared to the velocity fluctuations. Other examples of mea- 
surements of individual fluctuating flow quantities at supersonic and hyper- 
sonic speeds were reviewed in reference 27. It does appear that reasonably 
accurate measurements of the fluctuating, flow Variables in compressible flows 
can be achieved. 

C. Experiments Supporting Turbulence Modeling Development 
in Two-Dimensional Flows 

In order to illustrate how the experimental foundation for turbulence 
modeling might be established, it is helpful to focus on a particular problem. 
Consider the example of transonic flow about an airfoil as sketched on the top 
of figure 12. Three characteristic regions are identified that have certain 
distinct features as listed in the left-hand column. Actually, such features 
are typical of those found in a variety of aeronautical applications, so tur- 
bulence models developed to solve this particular problem will obviously apply 
to more general cases. Building-bloCk. experiments first selected to evaluate 
and guide the model development for each of the regions are sketched in the 
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n...xt column. The numbers In pnfenth«8ea below those sketches re for to the 
references which contain the experimental data. 

An attempt was made to select these experiments on the basis ul the 
requirements listed previously. However, in order to establish as wide a data 
base as possible* not all the experiments clearly meet the complete require- 
ments. As will be discussed in the next section, model development has pro- 
ceeded to the point where unresolved issues regarding the details of the avail- 


able experimental data have surfaced and the need for additional experiments 
to resolve these issues is clear. These issues and sketches of additional 


experiments complete the information in figure 12. 

For the case of fully turbulent attached flows over airfoils where the 
boundary -layer approximation to the full Navier-Stokes equations applies, 
pressure gradient and walL curvature encompass the important flow features. 
Obviously, if the Mach number was increased heat transfer would also be 
important. The first experiments chosen will mainly test the ability of the 
model to predict the effects of adverse pressure gradient. One of the-experi- 

ments also includes a favorable gradient ; 1 ^ 

In the second region, separation and reattachment may occur and the list 
of flow features is longer and obviously the problems are more complex. For 
example, the inviscidrv.iscid interaction affects the whole flow field and the 
complete Navier-Stokes equations may be required to solve the problems. The 
flows chosen for model development all have rather severe interaction and test 
the ability of the model to predict separation extent and reattachment char- 
acteristics. The trailing-edge flow region has fewer experiments that provide 
either verification or guidance for modeling. Those first chosen test the 
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ability of Che model to predict separntiftrt in the presence of adverse pressure 
gradients and how a model used in attached flows might U extended with wake 
or U railirig-edg6 flows. 

The additional experiments sketched in the last column might be proposed 
so as to obtain answers to the unresolved issues arising from a study of 
modeling for the building-block experiments as discussed in the next section. 


IV. STATUS OF TURBULENCE MODELING FOR COMPRESSIBLE FLOWS 


An experimental foundation for guiding and verifying turbulence model 
development was proposed in the previous section. The status of the various 
modeling concepts will now be reviewed by comparing computations that employ 
the models with the experimental data for the attached and separated flow 
cases . 


! 


A. Attached Flows 

1. Boundary-Layer Flows .with Zero or Mild Pressure Gradients 

Algebraic eddy viscosity or mixing length closure models have been suc- 
cessfully applied to compute compressible, fcero-pressure, and mi Id-pressure 
gradient boundary- layer flows. See, for example, chapter 9 of reference 9 
where a large number of computations employing the baseline model described in 
section II is compared with experiment. For purposes of this review three 


additional examples are given to illustrate the success of the model and to 
point out where additional work may be necessary. 

A comparison of the experimental data of Allen on an axisynunetric 


ogive-cylinder body and a computation employing the baseline model is shown in 
figure 13. The computations were made with a boundary-layer program originally 
developed for laminar flow< 54 ) but modified several years ago by Marvin and 
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Sheaf fef so that turbulent flows could also be computed. Comparisons of the 
computations with the shift friction and velocity-profile data, obtained by 
three techniques, indicate the hind of agreement that can be expected when 
employing the base lift* model at low Mach numbers on bodies at adiabatic wall 
temperatures. Calculations with and without transverse curvature (I VC) -Co 
account for a thick boundary layer show little difference. 

At higher Mach numbers, modeling changes that account for heat transfer 
and density fluctuations must be considered. For Mach numbers up to 10, modi- _ 
fications to the baseline model apparently are sufficient to account for these 
effects. Shang^S) made this conclusion or the basis of comparisons with skin 

friction and heat transfer data. He accounted for the triple correlation, 

p VV T , neglected in the baseline model in an approximate way by multiplying 
equations (19) and (20) by the term [1 + (<p’>/p)l< The magnitude of the 
density ratio was obtained by solving the y-momentum equation for the turbu- 
lent fluctuating pressure and introducing its magnitude into the equation of 
state, along with a correlation for temperature fluctuations. Also he, as 
well as Others before him, found that the effects of heat transfer could be 
adequately accounted for by the introduction of a turbulent Prandtl number, 
equation (24), which could be expressed as a constant, Pr t = 0.9, or a variable 
value depending upon the location in the boundary layer. However, the solu- 
tions were not too sensitive to the choice of the Prandtl number and the value 
of 0.9 was recommended. Above a Mach number of 10, the adequacy of the base- 
line model or its modified form has not been established. Lack of sufficiently 
accurate data above M = 10 contributes substantially to this conclusion. 

Figures 14 and 15, taken from reference 55, illustrate the success of 
employing a modified baseline model to predict skin friction and heat transfer 


at hypersonic Mach numbers. At Mach 7*4* including the effects of the density 
fluctuations has a sliial.1 influence and the unmodified baseline model predicts 
the data quite well. At Mach number 10.5, the effects of density fluctuations 
are larger aiid the modified model apparently yields a better prediction of the 
skin friction data. In both these examples the turbulent Prandtl number was 
held constant at 0.9 and the predictions of heat transfer are conservative. 
Introducing a variable Prandtl number would result in Somewhat better agree- 
ment but apparently there is further need for study to resolve the question of 
adequately predicting the heat transfer. Also in these examples, the initial 
rise in skin friction subsequent to transition is not predicted very Well, but 
this result can be attributed to the influence of low Reynolds number on the 
Clauser modeling constant, 0.0168. (See, e.g., chapter 6 of ref. 9.) 

Uniformly distributed mass transfer at the wall is another problem of 
interest receiving greater attention now because^of the possibility of its use 
to reduce drag on aerodynamic surfaces and thereby decrease aircraft fuel con~ 
sumption. Modification to the baseline model to account for mass transfer at 
the wall .is accomplished by altering the value of A + in equation (22). An 
appropriate variation of A + with mass transfer is found through analysis of 
experimental data. This is illustrated in figure 16 where experimentally 
determined values of A + are correlated with a mass transfer parameter B/c£, 
where B = P w v w / p e u e c f is value of skin friction. The data 

are adequately represented by an exponential curve. 

Predictions of the skin friction and neat transfer employing this modifi- 
cation to A + have been made with the Matvin-Sheaffer^ 54 ) bound ary- layer 

program. Results are shown in figure 17 where the skin friction and heat 

transfer-normalized by their untranspired values are plotted versus a 


tUmoiWinnlesR mass transfer parameter. 


The skin friction decrease with 


Increasing mass addition is predicted reasonably well, the effects of cam 
possibility , or Mach number, are small both theoretically and experimentally. 


Recently, similar results for skin friction to 


6.6 


were presented 


( r, r, ) 


The heat transfer decrease is also predicted well to M * 0,7. However, not 


shown are results at higher Mach number where the available data indicate a 


far smaller decrease in heat transfer than is predicted. This latter aspect 


is discussed in reference 63 and still requires further theoretical and 
experimental study. 


2. Boundary-Layer Flows with Severe Pressure Gradients 

The status of turbulence modeling for compressible boundary layers 
encountering rather Severe pressure gradients will be discussed by considering 
the attached flow experiments introduced in figure 12. These experiments and 
thf ir test conditions are summarized in figure 18 taken from reference 20, 

Attempts to predict the experimental data by using the baseline algebraic 
turbulence model introduced previously gave results that were not entirely 
satisfactory. However, in some instance^a reasonable prediction was 

achieved as will be shown subsequently* The poorest prediction was obtained 
for the compression ramp data of Sturek and Danberg^ 49 ^ as illustrated in 
figure 19, taken from reference 12, and serves as an example to show how 
modifications to the model were employed to improve the situation. In this 
'.nstance the baseline model predicts no rise in skin friction in contradiction 
to the data. Modifying the baseline model coefficients A + or A + and k 
(see eqs . (19) and (22)) as suggested for incompressible flows with pressure 

gradient^ 9 ’ 64 ’ 65 ) did not improve the prediction for all of the flows and in 
some cases produced poorer results. 


fha curve lab led X * 0, where X ** J,/<S, in figure 19 shows such a result 
when A + and k are modified to account for pressure gradient, However, 
Horstman was able to improve the predictions with the baseline model by_.in.tro-- 
ducing ah additional empirical relaxation parameter (X) in the modification 
for A + and k. For example, the curve labeled X = 7.5 shows excellent 
agreement with the data. Results. similar to this were obtained for each of 
the experiments in figure 18. However y the disadvantage of such an approach 
was that each experiment required a different relaxation scale to effect a 
good prediction and the available data were not sufficient to determine a 
universal correlation of relaxation scale length in terms, of. pressure gradient 
that would encompass the important parameters of Mach number and Reynolds 
number. 

In contrast to the former results, a recent study, has shown that the two- 
equation model of Wilcox^ 17 ) and the full Reynolds stress model of 
Donaldson^ 19 ) introduced in Section II provide overall a better means of pre- 
dicting these adverse pressure. gradient flows. Furthermore, no adjust- 
ments to modeling constants were required. 

To illustrate this finding, examples of the results taken from refer- 
ence 20 are shown in figures 20, 21, and 22 for the adverse pressure gradient 
flow of Peake et al.( 4 °) In figure 20, the mean flow pressure and velocity 
are shown along with the skin friction and form factor, H. The baseline model 
prediction is labeled Marvin-Shaef fer and the Wilcox and Donaldson .models are 
labeled Wilcox-Traci and ARAP, respectively. The .rise in skin friction is 
predicted very well by the tub-equation and Reynolds stress models. 

Typical. of the baseline model, the predicted increase lags the data 
because no additional production of turbulence Other than that introduced 
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through chnngos in the mean flow Is accounted for, .The axial distribution of 
the shape factor IS reproduced reasonably well, by all the models, but the 
baseline model, predicts a consistently higher magnitude. The velocity and 
shear stress profiles for two axial locations, one after a run of adverse 
pressure gradient and a second far downstream in the constant pressure region, 
arc shown compared with the predictions lft figures 21 and 22. Ihe velocity 
profiles are predicted very well by the Wilcox and Donaldson models. The 
shear stress profiles are predicted (overall) better by the multiequatron 

models . 

The advantage of the models employing additional differential equations 
to describe the behavior of the turbulence results from the expectation that 
they are capable of predicting a wide variety of flows without modification of 
constants and that they might provide the best means for calculating through 
separated regions such as those discussed in the next section. 

B. Separated Flows 

In the example proposed in figure 12 at the end of section III, separation 
and reattachment in the region of a shock wave interaction with a boundary 
layer were characteristic of the flow features that would have to be predicted 
by a successful turbulence model. Work directed toward model development for 
this problem has been under way for several years, both experimentally and 
computationally, and the progress on both aspects will be reviewed now...... 

1. The Building-Block Experiments 

The building-block experiments now being used by the group at Ames 
Research Center to study the behavior of turbulent boundary layers undergoing 
separation and reattachment after their interaction with an incident shock 
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wave are described. Each experiment, in wlUhh there IS a significant coupling 
of the viscous and itiviscid flow fields, has a companion computer simulation 
that uSeS the complete time-averaged, Navier-Stokes equations and requires a 
large, fast computer for a. solution. The experiments, in various degrees of 
completeness at this time, cover a wide range of Mach numbers and Reynolds 
numbers. 

The first two experiments involve transonic flows and particular interest 
is directed toward shock-boundary-layer interactions on wings. Figure 23(a) 
Shows the experimental arrangement for a verification experiment being 
performed on a wing that spans the test section of a high Reynolds number 
channel recently built at Ames Research Center. The. facility operates in a 
blow-down mode, and the free-stream Mach number can be adjusted before or 
during tests by a translating wedge that acts as a downstream choke. The 
upper and lower walls are contoured so that strong shocks would not extend to 
the walls and choke the. flow. A thick circular-arc wing (fig. 23(b)) was 
chosen to allow local airfoil Mach numbers to achieve values where shock- 
induced separation would occur, but without having the shock extend more than 
about two-thirds of the distance between the wing and the outer wall. Surface 
pressures for various. Mach .numbers and Reynolds numbers have been obtained as 
well as surface skin friction at specific Reynolds numbers. 66) Still to 
be obtained are mean velocity data and more skin-friction data. 

The wing experiment resulted in both shock-induced and trailing-edge- 
indueed separation, depending on the free-stream Mach number. Data were 
obtained to Reynolds numbers, based on chord. length, of 17*10 6 . At intermedi- 
ate Mach numbers, some unsteadiness in the. flow field Occurred. The results 
are discussed in detail in reference 45. ReynOlds-number effects for both the 
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trail I iig otlgp rtnd shock- Induced Separation were not significant for wliig-chord- 

based Reynolds numbers beyond 10*10 G . 

Some results obtained When shock-induced separation occurred are giviMi in 
figure 24. Oil-flow patterns (lower portion of the figure) illustrate the 
two-diniensionality of the flow and the detail in the region downstream of 
separation. The shadowgraph view near the interaction clearly illustrates 
Shock-induced separation. The pressure ratios downstream of the shock are 
below C*, indicating that the flow is slightly supersonic and suggesting the 
presence of an oblique Shock (also apparent in the shadowgraph) . 

Figure 25 presents the airfoil pressures at several Reynolds numbers and 
the skin-friction measurements at a single Reynolds number. The skin-friction 
measurements were obtained recently with surface-mounted wire gages specifi- 
cally developed for this experiment. ( 2G » The location of separation, 
determined from the oil-flow photograph, is shown on the abscissa of the skin- 
friction plot. This location also corresponds to the location .of the knee in 
the pressure curve downstream of the shock. The comparison with the computa- 
tions is discussed subsequently. 

Figure 26 shows the physical arrangement of another transonic flow 

(no) 

experiment being used for code verification and model concept development.' 
These tests are also being conducted in the Ames High Reynolds Number Channel. 
Supersonic flow was developed at the entrance of an axisymmetric test section 
and a normal shock wave was positioned at a fixed location by adjusting the 
location of a downstream shock generator. The relative distance between the 
shock wave and the downstream shock generator was always about 1 m. Experi- 
mental verification data were obtained for Reynolds numbers, based on distance 
along the wall to the location of the shock wave, between 9*10 6 and 290*10^. 
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Over this range rtf Reynolds number, Mach number varied between 1.35 and 1.45 
because of the differences in wall boundary- layer growth. With this arrange- 
ment, data are also being obtained at constant shock Math number by allowing 
the shock position to vary along the tube Surface when the Reynolds number is 
varied. A complete flow documentation, including turbulence measurements, 
has been performed at a Reynolds number of 37*10 6 where the corresponding Mach 

number ahead of the Shock wave was 1.44. 

Examples of the data are presented in figures 27 and 28. Additional 
results, including velocity profiles and velocity fluctuation data, are given 
in reference 29. The shock position is located at x/6 u = 0. The pressure 
rises rapidly downstream of the shock wave, causes separation, and then 
increases gradually. The corresponding Skin friction is reduced ahead, 
reverses sign in the separated region, and then increases thereafter. 

The shear Stress data (fig. 2.8) were obtained with the backed cross-wire 
anemometer described previously. The maximum shear stress in the boundary 
layer Shows a significant increase after the shock wave and then decreases 
downs cream toward the value expected for a mild adverse pressure gradient. 
Downstream of the shock wave beyond the effective boundary layer, Some mea- 
surable Shear was evident. It reverses sign because the mean velocity pro- 
files are retarded, probably resulting from a coalescence of compression waves 
ahead of the shock wave/ 2 *) In .this flow, the height of the separated zone 
was very small and no details were measured in the reversed-flow region. 

Documentation of the mean flow over a wide range of Reynolds numbers is 
now under way. These data will be used to assess the ability of the turbulence 
model to predict Reynolds-number effects on such factors as interaction 


lengths, pressure riSe to separation, and the effects of Math dumber on Intlp- 
iertt separation. 

One experiment being u§ed to verify the computations in the supersonic 
regime is the adiabatic shock“impinging flow originally reported by Reda and 
Murphy. (49) Figure 29 is a schematic of the test arrangement. Shock waves of 
varying strength were impinged on a Mach 3 boundary layer developed along a 
tunnel wall. For some tests, the Shocks were strong .enough to separate the 
flow on the tunnel Wall. Mean-flow profiles and surface-pressure data were 
reported originally. Since the original Work was completed, skin-friction 
data have been inferred from the profile data^ 67 ^ and shear distribution 

r ? 6 ) 

upstream and downstream of an interaction without separation were measured. 
Currently, measurements throughout a separated region, such as velocity fluc- 
tuations, shear distributions, and surface skin friction, are being documented. 

The surface pressure, normalized by the upstream total pressure, and 
indirectly inferred wall shear stress for a separated flow case are shown in 
figure 30. Separation. and reattachment points from oil-flow photographs are 
indicated. The pressure data have a plateau near the separation. Data for 
various shock strengths < 42 ) will be used to assess the ability of the Navier- 
Stokes codes to predict the onset of separation. 

At hypersonic speeds, the axisymmetfic shock expansion boundary-layer 
interaction flow described in reference 51 is being used to guide turbulence 
modeling concepts. The experiment is sketched in figure 31. The leading edge 
of the shock generator was varied between 5° and 20° . Measurements in the 
axial direction were made in finely spaced steps by traversing the shock 
generator in the axial direction during the tests. Complete flow documenta- 
tion, including turbulence measurements across the flow field, is available 


42 


for shock generator angles of 7.5°, where unseparated interactions occurred, 
and at 15°, whore fully separated interactions occurred. Surface pressure, 
skin friction, and heat transfer values ate available for other generator 
angles^ 42 ) and will be used to Verify the ability of the turbulence model to 
predict the effects of shock strength on separation at hypersonic speeds. 

The Surface pressure, Skin friction, and heat transfer for separated flow 
are shown in figure 32. Separation and reattachment points obtained from 
interpreting pitot-pressure measurements on probes near the surface facing 
both upstream and downstream are shown. The pressure increases through the 
interaction region, reaches a plateau near Separation, and then rises further 
after reattachment. The skin friction decreases, then rises downstream of 
reattachment; the heat transfer rises continually, the decay in pressure, 
skin friction., and heat transfer downstream of the interaction results from 
the expansion fan emanating from the comer of the generator. 

Fluctuating turbulence properties, for use as guides in modeling changes, 
were measured across the boundary layer at four locations through the inter- 
action. (^) an example of the measurements is shown in figure 33, where the 
shear stresses at the four measurement locations are plotted. It was not pos- 
sible to obtain shear measurements in .the_reversSd-flow region of the separa- 
tion bubble at the Second measurement Station so the expected trend has been 
sketched. The measurements show many of the same features as for the tran- 
sonic experiment Shown previously. Through the interaction, the maximum level 
of shear stress near the Separated region increases significantly. Downstream 
of reattachment, (x - x 0 )/6 Q > -1, the shear profiles do not differ in shape 
from those usually found in zero Or slightly favorable pressure gradients. 
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Turbulence memory apparently persists only for about five boundary^ Layer 
thicknesses beyond the initial rise in pressure.^ 5 ) 

The supersonic corner flow experiment ^ 52 ^ sketched in figure 34 is the 
Last of the shock interaction experiments to be discussed. Mean flow measure- 
ments at two ramp angles, 20° and 24°, have been documented and turbulence 
measurements are planned for the forthcoming year for the 20° ramp angle. 

Data have been obtained for an extremely large range of Reynolds number, based 
on boundary-layer thickness: 5x10° < Re^ < 7xl0 6 . Typical results showing 

surface pressure and skin-friction data for the 20° ramp angle taken from 
reference 52 are reproduced in figure 35, Features common to all the super- 
sonic interactions are evident; a pressure plateau and a rapid decrease in 
skin friction followed by a rise after separation, 

2, Companion Computations Using Various Turbulence Models 

Each of the experimentally determined flow fields described above have 
been programmed for numerical simulation on a CDG-7600 computer. Each program 
uses the mass-averaged forfn of the time-dependent Navier-Stokes equations and 
solves them with the MacOormack time-marching explicit scheme with split- 
ting^ 68 ^ and employs the more time-efficient algorithms proposed by 
MacCormack. ^ 69 ^ The complete Navier-Stokes equations were used because the 
viscous and inviscid flows interact significantly and it was-believed that use 
of approximate solution techniques would only introduce uncertainty when eval- 
uating the ability of the turbulence models to predict measured trends with 
Mach number and Reynolds-numb er variations. Each of the programs and solution 
techniques has been reported on separately and their details will not be 
addressed here. Solutions for each of the flows described above have been 
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made employing the baseline and relaxation models dismissed previously. Three 
of the flews have been predicted with the One-equation model and the others 
will be attempted in the near future. Eventually, it is likely that all the 
flows will require at least a two-equation model before reliable predictions 
can be expected. The results of the computations ate discussed next. 

The data from the transonic-wing experiment are compared in figure 24 
with the computations performed by Deiwert.^ 70 ^ With the baseline model, the 
overall features of the flow field are predicted reasonably well. Tor 
example, the trend of increasing airfoil peak Mach number with Reynolds number 
is predicted, as is separation due to shock interaction. Downstream of the 
shock, where the flow separates, the pressure recovery is overpredicted 
because the predicted shock wave is normal, whereas the experimental shock 
wave is oblique. A comparison of the predicted and measured skin friction at 
Re = 10 7 further illustrates the differences .in that region as the location 
of predicted separation is downstream of that location determined from the 
oil-flow photographs. Use of the relaxation concept (eq. (25)) with one 
boundary- layer thickness for L (relaxation length) tends to shift the loca- 
tion, of the shock. wave closer to the leading edge and reduces, the peak Mach 
number and resulting shock strength; nonetheless, the pressure recovery is 
still overpredicted. 

Solutions using other choices for L did not improve the results. For 
Mach numbers at which the flow is not separated but at which viscous effects 
are still important, the computations using the baseline model compare very 
favorably with the data, except. at the trailing edge where extensive separa- 
tion occurred. ( 7 °) Because solutions to this problem take between 30 and 
60 min of computation time, models employing differential equations for 
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turbulent kinetih energy or shear have not been attempted. These will be 
added in the near future, however, if other, more simple problems show 
improvement with employment of these models. 

Comparisons of the computations with the data from the normal shock-wave 
experiment are presented in figures 27 and 28. The pressure predictions 
employing all the models agree well with the. data. The Skin friction is 
underestimated downstream, using either the baseline or relaxation models. 
However, the one-equation model of Rubesin does significantly better with 
regard to predicting the downstream skin friction and is considered a major 
improvement . The latter results were reported recently by Viegas and 
Coakley. The shear distributions obtained using the. relaxation model are 
compared with the data in figure 28 at three locations downstream of the shock 
wave. Similar results were obtained with_...the baseline model. The computations 
fail to predict the significant increase in shear at the first station and, 
evidently, this continues to affect the prediction downstream. Better com- 
parisons are expected with the one-equation model, but they were not available 
at the time this paper was written. 

Comparisons for the oblique shock ihteraction using the baseline and 
relaxation turbulence models are shown in figure 30. For the oblique shock 
interaction, the baseline model predicts the overall pressure rise reasonably 
well and also predicts separation. However, no appreciable upstream influ- 
ence in the pressure and no plateau are indicated in the calculations because 
the size of the separation bubble is underpredicted. Introducing the relaxa- 
tion model with L = 5 boundary-layer thicknesses corrects this deficiency 
because it lowers the effective eddy viscosity near the separation, thus 
increasing the size of the separated zone, which in turn introduces a plateau 
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In the pressute-tise Curve. The location of separation is predicted with t lie* 
relaxation model, but the reattachmertt location is too far downstream. 

The comparison with the wall shear downstream of reattachmertt where the 
boundary layer thins may appear better than is the actual case because these 
computations were made with a relatively coatSe grid, and in the downstream 
locations the first computation point away from the wall was in the logarithmic 
region of the turbulent boundary layer. Since the method used to obtain 
derivatives at the wall requires calculated points within the sublayer region, 
the.. predicted, shear is expected to.be lower than the measurements, and the 
results would be similar to those for the transonic, normal shock comparisons* 

In figure 32, the predicted results and experimental data for the hyper- 
sonic shock interaction are compared. As With the supersonic oblique Shock 
interaction, the baseline model prediction reproduces the overall trends in 
the data but is deficient in the separated region. Recent calculations by 
T. J. Coakley^ 71 ) using the one^-equatipn model of Rubesin improves the situa- 
tion. The height of the separation was greater and the resulting upstream 
influence leads to. a plateau .in the pressure. However, Coakley had to alter 
the algebraic length Scale in the model to affect the changes. The .location 
wheie skin friction begins to decrease is predicted, but the extent of the 
separation in the axial direction is still overpredicted. The heat transfer 
is underpredicted throughout the interaction region with either model. 

The Supersonic corner. flow has been predicted using the baseline, relaxa- 
tion, and one-equation turbulence models. The latter computations were 
recently made by Horstman. fhe results- are compared with the experimental 
data in figure 35. The baseline model predictions show no plateau in the 
pressure distribution as a result of the small predicted separation zone, and 



the Skirt friction downstream Of rOat tachiiiertt is underpfedicted , The relnxa- 
tiort model does increase the size Of the separation zone and a plateau in 
pressure is predicted, but the skin friction downstream of reattachmertt is 
still underpredicted . The one-equation model significantly improves the pre- 
diction of skin friction and does reasonably well in the prediction of the 
overall pressures. Moreover, a comparison recently made by Horstman of the 
velocity profiles downstream of reattachment with the computations shows a 
marked improvement when the one-equation model is used. 

Prom an examination of all the comparisons, the following general obser- 
vations can be made. In most cases, the pressure rise can be predicted 
reasonably well with a simple eddy viscosity description for the turbulence. 

The exception is for the transonic Wing where the shock-induced separation was 
very large and extended from the foot of the shock beyond the trailing edge. 

The algebraic eddy viscosity concept is deficient in providing the proper 
details .within the separated zone and this, in turn, apparently affects the 
entire flow field when the Separation zones are large. 

Introducing relaxation tends to improve the situation somewhat because it 
decreases the eddy viscosity in the region of separation, resulting in an 
increase in the height of the separation. However, the next decrease in eddy 
viscosity persists beyond the separated zone even for short relaxation lengths, 
and t he codes underpredict the skin friction downstream. Interpretation of 
the modeling experiments indicated the need far an increase in effective eddy 
viscosity downstream of reattachment. An improved turbulence model that 
effectively accomplished this is the one-equation model of Rubesin. Further 
improvement to this model can be made by altering the algebraic length scale 
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as shown by Viegas and Coakley. < 71 > In fact, their results suggest that a 
twO-equation model will probably provide an even better prediction. Still 
needed, however, is a better model for the flow in the separated region where 
the flow is reversed; but the lack of data from such regions makes the guid- 
ance of model changes difficult. 

V. CONCLUDING REMARKS 

A complete Concept of turbulence modeling has been described. The ele- 
ments - experiments, modeling concepts, and computations - were each reviewed 
to provide an up-to-date status of the problem. Experiments were classified 
as verification or modeling types. Verification experiments were described as 
those that measured mean flow quantities and they would provide a means for 
assessing the ability to predict correctly over the flight ranges of Mach and 
Reynolds number. On the other hand, modeling experiments would require mea- 
surements of the turbulent flow itself in Sufficient detail to guide and 
verify the turbulence models actually used in the computation codes. 

Recent techniques used to measure turbulent properties in compressible 
flow were reviewed and the results are encouraging enough to permit the con- 
clusion that modeling experiments can be performed at various Mach numbers. A 
set of experiments to evaluate and guide turbulence model development for two- 
dimensional flows was introduced and comparisons of computations employing 
various turbulence models were made. The results of the comparisons showed 
that algebraic eddy viscosity models were sufficient to predict most attached 
flows with zero or mild pressure gradients up to Mach number 10. 

Models that employ additional differential equations to describe turbu- 
lent kinetic energy or shear stress are required in the Case of severe 
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pressure gradients . Flows with separation on ^attachment are also predicted 
better With models uSirtg additional differential equations, but the number of 
examples is limited at the present time and further comparison will have to be 
made before a predictive model is available for those flows. 
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FIGURE CAPTIONS 

Figure 1,- Process of physical modeling for computer code development. 

Figure. 2.- Status of turbulence modeling. 

Figure 3..- Mach and Reynolds number domain for aerospace vehicles. 

Figure 4.- Mach and Reynolds number domain for experiments compared with that 
for aerospace vehicles. 

Figure 5.* Comparison of estimated supersonic shear stress with incompressible 
measurements and the empirical prediction of Maise and McDonald. (Taken 
from ref. 27 . ) 

Figure 6.- Probes developed for high-temperature applications. (Taken from 
ref. 28.) 

Figure 7.- Supported crossed hot-wire probe. (Taken from, ref . . ) 

Figure 8.- Comparison of calculated and measured velocity-density sensitivity 

ratios. (Taken from ref. 31.) 

Figure 9.- Comparison of calculated density and velocity sensitivities and 
measured mass-flux sensitivity. (Taken from ref. 31.) 

Figure 10.- Shear-stress measurements in compressible flows. 

Figure . 11.- Normalized rms velocity and density fluctuations across a com- 
pressible turbulent boundary layer. 

Figure 12.- Model development for two-dimensional flows. 

Figure 13.- Comparison of baseline turbulence model computation and data of 

Allen. ( 53 ) 

(a) Skin friction. 

(b) Velocity profile, x - 2.08 ft. 

Figure 14.- Effect of density fluctuation of C f and St _at H q = 7.4. 

(Taken from ref. 55.) 
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Figure 15,“ Effect of density fluctuation on and St at a 10.5. 

(Taken from ref. 55.) 

Figure 16.- Effects of mass transfer on baseline model.. constant , A + . 

Figure 17.- Comparison between calculations using a modified baseline model 
and experiments for transpired boundary layers. 

Figure 18..- Experimental test flows used as standards of comparison for the 
turbulence models. (Taken from ref. 20.) 

Figure 19.- Comparison of measured and computed skin-friction distributions 
obtained using a turbulence lag model. ( 12 ' 

Figure 20.- Streamwise distributions of boundary-layer characteristics, experi- 
ment of Peake, Brakraann, and Romeskie. (Taken from ref. 20.) 

(a) Boundary- layer edge conditions. 

(b) Local skin friction coefficient, based on free-stream dynamic 
pressure. 

(c) Local form factor. 

Figure 21.- Velocity and turbulent shear stress profiles after run of adverse 

pressure gradient, data of Peake etal.s x = 52 cm. (Taken from ref. 20.) 

(a) Mean velocity profile. 

(b) Turbulent shear stress profile. 

Figure 22.- Mean velocity and turbulent shear stress profiles far downstream 
in Peake et al. experiment, x = 67.3 cm. (Taken from ref.. 20.) 

(a) Mean velocity profile. 

(b) Turbulent shear stress profile. 
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Figure 23. - Experimental arrangement for a transonic wing undergoing shoek- 

induced separation. 

(a) Overall arrangement. 

(b) Wing design. 

Figure 24.- Shock-induced boundary-layer separation experiment on a biconvex 
circular-arc airfoil; Re » 10.3x10®, M to “ 0.786. 

Figure 25.- Results from the. shock-induced separation experiment on a biconvex 

circular-arc airfoil. 

Figure 26.- Experimental arrangement for a normal shock-wave, turbulent 
boundary-layer experiment with separation. 

Figure 27.- Results from normal shock-wave experiment; Mr, 1*^4 and 
Re x = 37x10®.. 

Figure 28.- Shear-stress profiles from the normal shock-wave experiment, 

= 1...44 and Re x = 37x10®. 

Figure. -29.- Experimental arrangement for a supersonic oblique shock-wave, 
boundary-layer interaction experiment; M«, - 3.0 and Re^ * 9.73xl0 5 . 

Figure 30.- Results from the supersonic oblique. shock-wave, boundary-layer 
interaction experiment; M* * 3.0 and Re^ * 9.73xl0 5 . 

Figure 31.- Experimental arrangement for an axially symmetric hypersonic 

oblique shock-expansion, boundary-layer interaction experiment; M* * 6.9 
and Re x = 13x10^* 

Figure 32.- Results from the hypersonic oblique shock-expansion, boundary- ... 
layer interaction experiment; = 6.9 and Re x = 13x10®. 

Figure 33.- Shear-stress distributions for the hypersonic oblique shock- 
expansion, boundary-layer interaction experiment; = 6.9 and 
Re x = 13x10®. 
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Figure 34.- Sketch of 8- by 8-in. supersonic wind-tunnel and test model. 
(Taken.. from ref. 52.) 

Figure 35.- Results from the supersonic ramp flow experiment; 3, a 

and Re* * 
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Figure 2.- Status of turbulence modeling. 
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Figure 8.- Comparison of calculated and measured velocity-density sensitivity 

ratios . (Taken, from ref.,..__31 . ) 
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Figure 12.- Model development for two-dimensional flows 
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(c) Local form factor. 
Figure 20. - Concluded. 
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Fi cure 21.- Velocity and turbulent shear stress profiles after run of adverse 
pressure gradient, data of Peake et al.: x = 52 cm. (Taken fromx.cf, 20.) 
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(a) Mean velocity profile. 

22.- Mean velocity and turbulent shear stress profiles far downstream 
in Peake et al. experiment, x = 67.3 cm. (Taken ftom ref. 20.) 
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DIMENSIONLESS CHORD 
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Results from the shock-induced separation experiment on 
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Results from normal shock-wave experiment; 
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Figure 31.- Experimental arrangement for an axially symmetric hypersonic 
oblique shock-expansion, boundary- layer interaction experiment; M, - b.9 
and Re x = 13*10 6 . | 













